Semileptonic Decays of B** mesons into D s \ and D* s2 in HQET 

Long-Fei GaiJ*] and Ming-Qiu Huang 

Department of Physics, College of Science, 
National University of Defense Technology, 
Changsha, Hunan 410073, People's Republic of China 
(Dated: December 23, 2011) 

In the framework of the heavy quark effective theory, the leading order Isgur-Wise 
| functions relevant to semileptonic decays of the orbitally P-wave excited B s meson 

states B** , including the newly found narrow £? s i(5830) and -B* 2 (5840) states, into 
the ( I? s i(2536), L>* 2 (2573) ) doublet are calculated from QCD sum rules. With these 



universal form factors, the decay rates and branching ratios are also estimated. 
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O . I. INTRODUCTION 

The investigation of the semileptonic b — > c processes is an important source for the de- 
termination of the parameters of the standard model, such as Cabibbo-Kobayashi-Maskawa 
matrix element \V c b\. They also provide valuable insight in quark dynamics in the nonper- 
^ . turbative domain of QCD. So a considerable amount of work has been carried out in the 
£N) ■ semileptonic decays of ground-state B and B s mesons both experimentally and theoreti- 
cally during the last two decades With the running energies of the colliders enhanced 
and the precision of the detectors improved, more and more excited bs mesons have been 
observed. In 2008, two orbitally excited narrow bs mesons, P s i(5830) and _B* 2 (5840), were 



observed by the CDF and DO Collaborations 10, This has inspired a lot of interest in 

5-H ' 



their various decay modes 12 
our consideration. 



151 ] . among which the semileptonic decay modes are under 



The heavy quark effective theory (HQET) 16] has been proved to be an useful theoret- 



ical tool to describe the spectroscopy and matrix elements of mesons containing a heavy 
quark. In the heavy quark ttiq — > oo limit, the spin and parity of the heavy quark decouple 
from those of the light part. So the heavy-light mesons can be categorized according to 
the total angular momentum of the light degree of freedom j q . Combining j q with the spin 
of the heavy quark yields a doublet of heavy-light meson states j = j q ± 1/2. For S'-wave 
ground states with j q = 1/2 (L = 0), two states of negative parity form the H doublet 
(0~, 1~). For P-wave excited states with j q = 1/2 or j q = 3/2 (L — 1), four states of posi- 
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tive parity form the S doublet (0 + , 1 + ) and the T doublet (1 + , 2 + ). The four P-wave states 
of bs system are usually denoted as B** and the newly observed -B s i(5830) and £>* 2 (5840) 
mesons are considered as members of the T doublet. 

The QCD sum rule method |l7| has been proved to be a successful nonperturbative 
approach which was widely used to investigate the properties of various types of hadrons 
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18H2JJ. In our previous papers, we have employed it to investigate the semileptonic de- 
cays of ground-state B mesons into some highly excited charmed meson doublets in the 
framework of HQET [9J, and then extended the discussion to the case of the semileptonic 
decays of S and T doublets of bs system into H and S doublets of the cs system [15]. 
We find that the semileptonic decay widths of a doublet of bs system into the same dou- 
blet in the cs system are significantly enhanced rather than those of its transitions into 
a different doublet. So it is reasonable to expect that the semileptonic decay widths of 
(5,1(5830), 5*2(5840)) into (D al (2536), L>* 2 (2573)) should be comparable with the previ- 
ous decay modes. In order to verify this point, we investigate these semileptonic processes 
in this paper. As a byproduct, we also investigate the semileptonic decays of the other B** 
doublet into the same final states. 

The remainder of this paper is organized as follows. After an introduction, we derive 
the formulas of the weak current matrix elements in HQET in Sec. [Ill Then we deduce 
the three-point sum rules for the relevant universal form factors in Sec. IIHI In Sec. [TV] we 
give the numerical results and discussions. The decay rates and branching ratios are also 
estimated in the final section. 



II. PARAMETRIZATION OF DECAY MATRIX ELEMENTS IN HQET 

The semileptonic decay rate of a B s meson transition into a D s meson is determined by 
the corresponding matrix elements of the weak vector and axial- vector currents {V^ = cy^b 
and N 1 = 07^756) between them. These hadronic matrix elements can be parametrized in 
terms of form factors. In HQET, the classification of these form factors has been simplified 
greatly. For the transitions between the S doublet of the bs mesons and the T doublet 
of cs mesons, the weak matrix elements can be parametrized in terms of one Isgur-Wise 
function at the leading order of the heavy quark expansion. But for the transitions between 
(.8,1 (5830), B; 2 (5840)) and (L> sl (2536), L>* 2 (2573)), we need two Isgur-Wise functions to 
parameterize the corresponding hadronic matrix elements at the leading order of the heavy 
quark expansion. 



According to the formalism given in Ref. [22|, the heavy-light meson doublets can be 
expressed as effective operators. For the processes we considered, two heavy-light meson 
doublets S and T are involved. The operators Pq and P[ that annihilate members of the 
S doublet with four-velocity v are, in the form, 

S v = l -^[P[^ + Po\- (1) 
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The fields P v and P*^ v that annihilate members of the T doublet with four-velocity v are 
in the representation 

= ^{P*^l» - \fl pU ^W - \lu{-f ~ «")]}, (2) 

where f = v ■ 7. At the leading order of the heavy quark expansion, the hadronic matrix 
elements of the weak currents between states in the doublets S v and T v > can be calculated 
from 

itfrh® = ^y^K^f rsf }, (3) 

while the corresponding matrix elements between states annihilated by fields in T v and T v > 
are derived from 

h$rh& = -Tr{(Uy)g*« - UvW^T^rT^}, (4) 

where h^, are the heavy quark fields in HQET, and T v > = 70^,70. v is the velocity of the 
initial meson and v' is the velocity of the final meson in each process. The Isgur-Wise form 
factors rj(y), and £2(2/) are universal functions of the product of velocities y{= v • v'). 

Here we should notice that each side of Eqs. (J3J) and (J4]) is understood to be inserted 
between the corresponding initial B s and final D s states. The hadronic matrix elements of 
B s0 (B' sl ) D sl (D* 2 )£j7 can be derived directly from the trace formalism (jHJ) and are given 
as 

(D„v, c w-Ar \B M ) = lu lW _ _ (y + 2)l/v + 3 „ v 

^/m Ba0 m Dal V 

+ i(y-l)v p v' T e^], (5) 



(^(^,6)1(^-^1^0^)) 



^OK^Ki - v)g^ + - ^ p <^ r ], (6) 



+ (y - 1)K + v'»)g^ + z(y + 1)^> T - <) 



^v p v' T e^ pT - iv^Xe^] , (7) 



-ie^> T -<)]. (8) 

The hadronic matrix elements of — )■ D s i(D* 2 )iV are calculated similarly from Eq. (@| 
as follows: 

y/m Bal m Dgl b 

+ (2y + l)u V - 3vV(t/" + ^)] 
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+ Uv)\-^y 2 - 1)9^ + (i - v W" + «") 

- 2{y 2 - l)v p g pp + 3(y - 2)«V(i/" + ^)]}, (9) 

^m Bal m Dal 6 

+ v'e^X -(y + 2) £ ^>; + u<r )] 

- £ 2 (y)[(y 2 - + + (y + 1)*/ V' Q ^X 

- 2(y + 1)K - v'")^^ ~(y+ l)v p e^v a v' a }}, (10) 

y/m BBl m D . 2 V6 

+ g^ p e PaaT v a v' T + g pa e ppaT v a v' T + v'^e^Va - v' t e f,p< " x v' a 

- v^e^ paa v a + v p e^ aa v' a + v a e^ pa v a + (y - l)e^ pa ] 
+ &(y)[v a e^ pa v' a - v a v ,l3 e ppaT v a v' T 

- 2v r7 v' p e^ paT v a v' T - yv c e^ pa v' a - (y - l)v° 'e^ pa v a }}, (11) 
(D: 2 (v> , e')\A p \B sl ( V 2 e)l l^ ^^^f, + 2 ^ (( + ^ _ ^ 

y/m Bsl m Dh V6 

- + + (y + ^)g pp g^] 

+ £ 2 (y)[-2(l + y)v a g pp v'P - (y + l)v a g pp v' p 

- (y + l)v p v a gP p + {y + l)v p v a g pp + 3v p v a v'^v' p ]}. (12) 

For the decays B* 2 — > _D sl (_D* 2 )£z7, the corresponding hadronic matrix elements are as 
follows: 

(D^v ^WlBU v ^ = _ ^ eap{Uy)[g ^r Va< + 2g « p e^ VaV > T + g^ £ « p - VaV > 
y/m B * 2 m Dsl V6 



+ v% a ^X + v' p e afipT v T - v p e afipT v' T + (y - l)e aPpp ] 

+ Uy)[{y - l)v ,a e^ pT v T + (y- l)v' a e^v' T 

- 2v p v' a e ppaT v a v' T - v p v' a e^ aT v a v' T }}, (13) 

(D sl {v', e')\A p \B* s2 {v,e)) = 1 {^ y )[ g ^ v '^ _ /VV + 2 ^((y + l)g^ - v p v ,p ) 

^m B:2 m Dsl ^6 

+ g a »((y + 1)/' - vPv'?)] 

+ 6(y)h(l + y^W' 3 - (y + l)g? p v' a (v p - v ,p ) 

+ 2(y + l)^Vt/° + St/WV^]}, (14) 



(D; 2 K,6Q|^|i?; 2 (t;,6)) 



^ p e a( s{Uy)9 a V P ^ - + v p ) + v p g^\ 

- &{y)v a v' a [g pp v'P - gP p {v' p + v p ) + v p g^}}, (15) 



5 



(DUv ^^lBU v ^ { ( } + } _ 6( ^V^« + v T )} 

y/m Bh m D * 2 

(16) 

In these matrix elements, e a (e' a ) is the polarization vector of the initial (final) vector 
meson and e a p (e' ) is the polarization tensor of the initial (final) tensor meson. The only 
unknown factors in the matrix elements above are r)(y), £i(y) and £2(2/) which should be 
determined by nonperturbative methods. In the following section, we will apply the QCD 
sum rule approach to estimate them. 



III. FORM FACTORS FROM HQET SUM RULES 



In order to apply QCD sum rules to study these heavy mesons, we must choose ap- 
propriate interpolating currents to represent these states. Here we adopt the interpolating 
currents proposed in Ref. |23| based on the study of Bethe-Salpeter equations for heavy 
mesons in HQET. Following the remarks given in Ref. jsj, we take the interpolating cur- 
rents that create heavy mesons in the T doublet and S doublet as 



l,+,3/2 



rfa/3 
2,+,3/2 



jt 

0,+,l/2 



-I, 



.-%) 



V2 



T a ^h v ^D tu s, 



-j=h v (-i)ft t s, 



J 



fa 



,1/2 



-i)p t s, 



(17) 
(18) 
(19) 
(20) 



where = D a — v a (v ■ D) is the transverse component of the covariant derivative with 
respect to the velocity of the meson. The tensor T al3 >^ v is used to symmetrize the indices 
and is given by 

T a ^ = \{gTi v + 97 i») - \gfgr, (21) 

where = g al3 — v a v 13 is the transverse part of the metric tensor relative to the velocity 
of the heavy meson. 

These currents have non-vanishing projections only to the corresponding states of the 
HQET in the iriq — > oo limit, without mixing with states of the same quantum number 
but different S[. Thus we can define one-particle-current couplings as follows: 



(H sl (v,e)\J« +t3/2 
(H. 



|0> 



s2\ u i c VI' v 2,+,3/2l w / — /2,+,3/2a/^* 2 ^ > ' 

(H sl (v,e)\J^ + ^ 1/2 \0) = fo,+,i/2^/m Ha0 , 
(H* s2 (v,e)\JZ +A/2 \0) = h t+ ,xf2yjm^e* 



for J p 
for J p 
for J p 
for J p 



r 

2" 

r 



(22) 
(23) 
(24) 
(25) 
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The decay constants /i,+,3/2, /2,+,3/2, fo,+,i/z, and /i,+,i/2 are low-energy parameters which 
are determined by the dynamics of the light degree of freedom. 

With these currents, we can now estimate the Isgur-Wise functions £i(y), £2(2/), and 
r)(y) from QCD sum rules. Let us consider £1(2/) and £2(2/) first. The jumping-off point is 
the following three-point correlation function: 

E^(u,u\y)=i 2 J rf^e 4 <*^^ 

= Si(w, w', 2/)£^ A) + H 2 (u;, u/, y)C a £* VA)y (26) 

where ^ = /i(t>')7 M /i(t>) and J^ 1 '' 1 ' ^ = /i(t>')7 M 7 5 /i(f ) are the weak currents, ^^3/2 is 

the interpolating current defined in Eq. ( ITT)) . Ei{uj,uj' ,y) (i = 1,2) are analytic functions 
in u = 2v ■ k and u' = 2v' ■ k', and are not continual when oj and u' locate on the positive 
real axis. k{= P — m^v) and k' (= P' — m c v') are the residual momenta of the initial and 
final meson states, respectively. The scalar functions Ei(uj, u , y) (i = 1,2) also depend on 
the velocity transfer y — v ■ v'. C^^^ and C^tyA) are Lorentz structures. 

To calculate the phenomeno logical or the physical part of the correlator (1261) . we insert 
two complete sets of intermediate states with the same quantum number as the current 
Ji +> 3/ 2 (a;) and isolate the contribution from the double pole at uj = 2A, uj 1 = 2A: 

s*">, co', y) = -= %f — ^—XMqr + Uy)qr\ + • • • > (27) 

(2A — u — zeJ(2A — uj ' — ze) s 5 

where "• • • " denotes contributions from higher resonances and continuum states, /i,+,3/2 is 
the decay constant defined in Eq. (|22|) . As we can see from the equations (1261) and (1271) . 
the pole contribution to Hi (a;, uj' , y) is proportional to the universal function £i(y) while 
that to E 2 (uj, uj\ y) is proportional to £ 2 (y). By isolating the different Lorentz structures, 
the sum rules for £i(y) and £2(2/) can be constructed from E\(u, u , y) and E 2 (uj,u\y), 
respectively. 

The theoretical side of the correlator is calculated by means of the operator product 
expansion. The perturbative part are usually expressed as a double dispersion integral 
in v and v plus possible subtraction terms. Therefore the theoretical expressions for the 
correlation functions in ( 126]) are of the form 

/pert / ' \ 
dudu- Pl ' 2 \' t ! V ~ 7 + subtractions + Hf^fw, u,y). (28) 
[U — UJ — 1£){U — UJ — IE) ' 

The perturbative spectral densities can be calculated straightforward from HQET Feynman 
rules. We consider only the leading order of perturbation here and the perturbative spectral 
densities of the two sum rules for £i(y) and £2(2/) are, 

C>> ^ y) = - 2^ (y+ ^4 -i^ K^ + "0 - 16m ^y + - + ^) 
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x ®{v)Q{v')Q{2yvv' - v 2 - z/ 2 ), (29) 

and 

"'•!/)=- ^ ( , + 1)7 4_ 1)5/2 {^ + - 1)(- 2 + O - 2(3, 2 - „ 

+ l)zV] - 64m s (y + l)[3yz/ 2 - 2vv'{2y 2 + 1) + 2yu' 2 }} 

x ®(u)Q(u')Q(2yvv' - v 2 - z/' 2 ), (30) 

respectively. Assuming quark-hadron duality, the contributions from higher resonances 
are usually approximated by the integrations of the perturbative spectral densities above 
some threshold. Equating the phenomenological and theoretical representations, the con- 
tributions of higher resonances in the phenomenological expression ( 1271) can be eliminated. 
Following the arguments in Refs. ll. UA\. we can not directly assume local duality between 
the perturbative and the hadronic spectral densities, but first integrate the spectral densi- 
ties over the "off-diagonal" variable z/_ = u — u, keeping the "diagonal" variable v + = ^p- 
fixed. Then the quark-hadron duality is assumed for the integrations of the spectral den- 
sities in u + . The integration region is restricted by the G functions above in terms of the 
variables v_ and u + and usually the triangular region defined by the bounds: < u + < u c , 
— 2y // |^"z/ + < V- < 2y/|^iz/ + is chosen. A double Borel transformation in u and u' is 
performed on both sides of the sum rules, in which for simplicity we take the Borel param- 
eters equal [l|, 0, 0] : T\ = T 2 = 2T. It eliminates the subtraction terms in the dispersion 
integral ( l28|) and improves the convergence of the operator product expansion series. Our 
calculations are confined at the leading order of perturbation. Among the operators in the 
operator product expansion series, only those with dimension D < 5 are included. For 
the condensates of higher dimension (D > 5), their values are negligibly small and their 
contributions are suppressed by the double Borel transformation. So they can be safely 
omitted. Finally, we obtain the sum rules for the form factors £i(y) and £2(2/) as follows: 

6(y)A 2 + ,3/ 2 e- {2Al ^ /2)/T ~ £ du + e-**[ui + m s (y + 1)^] 

{g 8 sa ■ Gs) + T 3y + 1 q, 



12 v 4T y 32(y + l) 2 n 

1 1 /"We 

6(y)A 2 + , 3/2 e- (2Al '+^ )/T =^2(^1)4 £ dv + e-^[34 + 2m s (y + iy + ] 

--7^ ± 4(-GG), (32) 

where {g s sa ■ Gs) = m^ss) with m 2 , = 0.8GeV 2 . 

The derivation of the sum rule for the Isgur-Wise function rj(y) is completely similar. 
Only now the correlation function we need to consider is 



1 



■• 2 ' d'xd'ze^ ^•HQ\T[J? +fl Jx)J$ v ) (0)jj +a/2 (z)|0> = S( W)W ',y)/%, (33) 
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where Jy^° ^ are also the weak currents, J" + 3 , 2 and Jo,+,i/2 are the interpolating current 
defined in Eqs. ( 1TT1) and (1T91 . Here only one Lorentz structure Cy A appears. By repeating 
the procedure above, we reach the sum rule for 7](y) as below: 

>f(i/)/o, + ,i/2/i. + , 3 /2e- (A "'"' +A ''«" )/T = 16 1 ?r2 ( _^ 1)3 J " *+e-*K + 2m,(y + 

+s^(v + iyj' + <*^(i-^) 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

Now let us evaluate the sum rules derived in the previous section numerically. First, 
we specify the input parameters into our calculation. For the vacuum condensation pa- 
rameters, we adopt the standard values: (qq) = — (0.24) 3 GeV 3 , (a s GG) = 0.04 GeV 4 , and 
(ss) = (0.8 ± 0.2) (qq). The mass of the strange quark is m s = 150 MeV. For masses of the 
initial B sl , B* 2 , B^ and B' sl mesons, we use M Bsl = 5829.4 MeV, M B * = 5839.7 MeV, 



M Ba0 = 5718 MeV [25J], and M B > si = 5765 MeV. For masses of the final D sX and D* 2 mesons, 
we use M Dsl = 2535.4 MeV and M D * 2 = 2572.6 MeV 

In order to obtain information of Isgur-Wise function £1,2(2/) and r](y) with less sys- 
tematic uncertainty, we can divide the three-point sum rules ( l3Tj) . ( 132|) . and (134"]) with 
the square roots of relevant two-point sum rules for the decay constants, as many authors 
did 0,H,B. This can-not only reduce the number of input parameters but also improve 



stabilities of the three-point sum rules. The two-point QCD sum rule we need here are 27] 
/i 2 +,3/ 2 e- 2A +^ /T ~ £ due-HS + 2my) - ^(ss) - ^GG)T. (35) 



and 

3 r^co y 
/o, + ,i/2e" 2A °' + ' 1/2/T =^ ^ dve~^(v A + 2m s u 3 - 6m> 2 - 12m 3 s u) - —m 2 (ss) 

x (1 - ^ + -^f) + -m 2 s (ss) - j^(-^GG). (36) 

After the divisions have been done, the Isgur-Wise functions £1(2/), £2(2/), an d r](y) depend 
only on the Borel parameter T and the continuum thresholds. The determination of the 
Borel parameter is an important step of sum rules. After a careful analysis, we find the sum 
rules for £i(y) and fe(y) have a common sum rule "window": 0.55 GeV < T < 0.65 GeV, 



which overlaps with that of the two-point sum rule fl35|) [27] . For the sum rule of r](y), 
we choose the "window" as 0.5 GeV < T < 0.65 GeV. Note that the Borel parameter 
in the three-point sum rules is twice that in the two-point sum rules. In the evaluation, 
we have taken 2.9 GeV < uj c < 3.1 GeV and 2.9 GeV < u c0 < 3.1 GeV. The regions of 



9 



>; 1 - 



0.95 - 
0.9 - 

0.85 - ( a ) 
0.8 L 



0.5 



a> c =2.9 GeV 
- co c =3.0 GeV 
-co =3.1 GeV 



T (GeV) 




FIG. 1: (a) Dependence o/£i(y) on Borel parameter T aty = 1. (b) Prediction for the Isgur-Wise 
functions £i(y) atT = 0.6 GeV. 




0.5 0.55 0.6 0.65 0.7 ' 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 



T (GeV) y 

FIG. 2: (a) Dependence of^iu) on Borel parameter T aty = 1. (b) Prediction for the Isgur-Wise 
functions £2(2/) atT = 0.6 GeV. 



;hese continuum thresholds are fixed by analyzing the corresponding two-point sum rules 
27]. Following the discussions in Refs. [l[ 24], the upper limit co c for u + in (I3T]) and 



(|32|) is just the same as that for v in (1351) . For w cl in Eq. (134"]) . it should be in the 
region + 1) — a/?/ 2 — l]w c o ^ ^ci ^ |(w c o + w c ). So it can also be fixed in the region 
2.9 GeV < U& < 3.1 GeV. The results are shown in Fig. [U Fig. [21 and Fig. EJ where we 
have fixed u c = 3.0 GeV in the two-point sum rule fl35l) and u c o = 3.0 GeV in Eq. fl36l) . 

The curves for £i(y), £2(2/), and 77(7/) shown in the figures above can be parametrized by 
the linear approximations 

£i(y) = £i(l)[l-4fe-l)], ^(1) = 1.03 ±0.03, p| = 0.75 ± 0.04; (37) 
Uv) = - pUv - i)], 6(1) = i-oo ± 0.03, p\ = 1.4 ± 0.1; (38) 
r){y) = r?(l)[l - pj(y - 1)], r/(l) = 0.45 ± 0.03, p% = 1.02 ± 0.03. (39) 



The errors mainly come from the uncertainty due to w c 's and T. It is difficult to estimate 
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(b) 



T (GeV) 



1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.< 

y 



FIG. 3: (a) Dependence ofr](y) on Borel parameter T at y = 1. (b) Prediction for the Isgur-Wise 
functions rj(y) at T = 0.58 GeV. 



the systematic errors which are brought in by the quark-hadron duality. Using the linear 
approximates for the universal form factors above, we can calculate the semileptonic decay 
rates of processes B s i(B* 2 ) — > D s i(D* s2 )lV and B s0 (B' sl ) — > D s i(D* s2 )tv. For this purpose, 
we have to derive firstly the formulas for the differential decay rates of these processes in 
terms of the Isgur-Wise functions £,i(y), £2(2/), and rj(y) from the matrix elements (|5l)-(|T6l) 
given in Sec. [Ill After some derivation, the formulas of the differential decay rates of the 
processes B s i(B* 2 ) — > D s i(D* 2 )£v appear as 

^(B sl -+ D A tu) = ^ l ^ l ^f;f L te(y) 2 [(r 1 2 + l)(26y 3 + 18y 2 + 21y - 11) - 2 n 

x (16y 4 + % 3 - I2y 2 -y + 43)] + 2&(y) 2 (y 2 - l) 2 [(r 2 + l)(13y - 4) 

- 2 n (% 2 - Ay + 5)] - 2Uy)Uv){y 2 - i)[(^ 2 + i)(26y 2 + h y - 13) 

-2r 1 (16y 3 -3y + 5)]}, (40) 
dY G 2 F \V c b\ 2 m 2 R rn 3 n * 

— (B sl D* s2 £v) = 34 x2 %3 tf 1 ^) ^ + + 54?/ + 159y + U ) " 2r2 

x (32y 4 + AOy 3 + 15Qy 2 + 25y + 17)] + 2&(y) 2 (y 2 - l) 2 [(r 2 + l)(23y 
+ 4) - 2r 2 (16y 2 + Ay + 7)] - 2^(y)^(y)(y 2 - l)[(r 2 + l)(46y 2 + 31y 
+ 13) - 2r 2 (32y 3 + 2Ay 2 + 27y + 7)]}, (41) 

dY G 2 j?\V cb \ 2 ml,* m 3 n 

-(B* s2 -+ D S1 £V) = 5 x 33 x [(r3 + 1)(4% + 5% + 15% + H) " 2rs 

x (32y 4 + 40y 3 + 156y 2 + 25y + 17)] + 2&(y) 2 (y 2 - l) 2 [(r 2 + l)(23y 
+ 4) - 2r 3 (16y 2 + Ay + 7)] - 2£ 1 ( 3 /)6( 2 /)(y 2 - l)[(r 2 + l)(46y 2 + Sly 
+ 13) - 2r 3 (32y 3 + 2Ay 2 + 27y + 7)]}, (42) 

dY G%\ V cb \ 2 m 2 R * m 3 n * 

-(B; 2 -> 7J S VT7) = 5 x 3 3 x ^ 57r 3 [("4 + 1)(7V + 66y 2 + 141y - 11) - 2r 4 
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x (4% 4 + 4(h/ 3 + My 2 + Ihy + 83)] + 2£ 2 (</) 2 (y 2 - l) 2 [(r 4 2 + l)(37y 

- 4) + r 4 (-48?/ 2 + 8y - 26)] - 2£ 1 (y)£ 2 (y)(y 2 - l)[(r 4 2 + 1)(74 2/ 2 

+ 29y - 13) - r 4 (96y 3 + 32y 2 + 26y + 26)]}, (43) 

while for the processes B s0 (B' sl ) — > D s i(D* 2 )£T/, they can be found to be 

^(B s0 -+ D sl iv) = G2AV f^f D - \v{y)\\y - i) 5/2 (y + i) 3/2 [(r 5 2 + l)(2y + 1) 

-2r 5 (y 2 + y + l)}, (44) 

^(b s0 -+ D a2 tv) = G2AV f^f DM \v{y)\\y - i) 5/2 (v + i) 3/2 K + my - 1) 

-2r 6 (3y 2 -y + l)], (45) 
-2r 7 (5y 2 -y + 2)\, (46) 

^(B' sl -+ D s2 iu) = G2AV ff^ Da2 Hy)\ 2 {y - if /2 {y + i) 3/ % 8 2 + i)(iiy + 1) 

-2r 8 (7y 2 + y + A)}, (47) 



where r$ (z = 1, • • • , 8) is the ratio between the mass of the final D s meson and that 
of the initial B s meson in each process, e.g., r\ = A/ ^ sl . The maximal values of y for 
these semileptonic processes are given in Table [IJ In addition, we need \V c f,\ = 0.04 and 

TABLE I: The maximal value of y for each process: y m ax = (1 + rf)/2n (i = 1, 2, • • • , 8). 





B s i 


B S 2 


B s o 


B' sl 


D sl eu 


1.36707 


1.36872 


1.34934 


1.3568 


D s2 eu 


1.35364 


1.35525 


1.33628 


1.34358 



Gf = 1.166 x 10~ 5 GeV -2 . By integrating the differential decay rates over the kinematic 
region 1.0 < y < y max , we get the decay widths of these semileptonic decay modes which 
are listed in Table [Til Although the widths of B s] and B* 2 have not yet been measured 
experimentally, they were estimated early in Ref. [28] to be around 1 MeV. Theoretically, 



their strong decays were investigated in Ref. [13]. As we know, the main decay modes 
of these excited B s mesons are strong decays. Therefore we can approximately take the 
strong decay widths as the total widths for an estimation of order of the branching ratios 
of these processes. In fact, the two-body strong decay widths of B s \ and B* 2 are computed 



to be 98keV and 5MeV in Ref. [131 ] . The masses of B s0 and B' sl are considered to lie 



below the thresholds of B*K and BK, so their main decay modes are isospin violating 
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TABLE II: Predictions for the decay widths and branching ratios 



Decay mode Decay width (GeV) Branching ratio 





-»• D sl £V 


(2.1 ±0.2) 


x 10 - 


14 


r\J 


10" 


10 


B s i ~ 


-> D* s2 iu 


(2.8 ±0.2) 


x 10 _ 


14 




io- 


10 


B s2 ~ 




(1.8 ±0.1) 


x 10- 


14 


r^j 


io- 


ii 


B* s2 - 


-> ^ 


(2.9 ±0.2) 


x 10- 


14 




10" 


ii 


B.so - 


->■ D sl £V 


(1.1 ±0.2) 


x 10" 


16 


r^j 


10" 


12 


B s o - 


+ D* 2 tv 


(1.0 ±0.2) 


x 10- 


16 




io- 


12 


B' sl ~ 




(0.8 ±0.1) 


x 10- 


16 


r\J 


io- 


12 


B' sl ~ 




(1.5 ±0.3) 


x 10- 


16 




io- 
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decays and radiative decays and are supposed to have width of about 100 keV [25[]. Using 
these widths, we estimate the orders of the branching ratios of the semileptonic decays (see 
Table ILT1) . Note that the semileptonic decay branching ratios of B s0 and B' sl should be 3 
orders of magnitude lower if their masses lie above the thresholds of B*K and BK. As 
we can see from Table HT| the semileptonic decay widths of (B sl , B* 2 ) into (D s i, D* 2 ) are 
significantly larger than those of (B s0 , B' sl ) into (D s i, D* 2 ), as was expected. The present 
precision of the experimental measurement of the branching ratio of the B s mesons has 
reached up to lO -7 ~ lO -8 [26]. Therefore, all the branching ratios of these processes are 
so small that it is difficult to find them in experiments. However, we can still expect that 
with the precision of the detector improved, the decays B sl — > D sl (D s2 )£v might be seen 
in the LHCb experiment. 

In summary, we have performed a study of the semileptonic decays of the orbitally P- 
wave excited B s meson states B**, including the newly observed -B s i(5830) and £>* 2 (5840) 
mesons into the T doublet of D** mesons, -D s i(2536) and _D* 2 (2573), within the framework 
of HQET. We employ QCD sum rules to estimate the leading-order universal form fac- 
tors describing these weak transitions. Different from the semileptonic processes of lower 
Heavy-light meson states, two universal Isgur-Wise functions £i(y) and £2(2/) are need to 
parameterize the hadronic matrix elements in the weak transitions B s i(B* 2 ) — > D s i(D* 2 )iV 
at the leading order of the heavy quark expansion. The predicted branching ratios of these 
processes are prohibitively tiny, so it is difficult to find them in experiments. The decay 
widths of B s i — > D s i(D s2 )£v are comparatively large so that we expect they might be seen 
in the future LHCb experiment. It is worth noting that the Isgur-Wise functions £i(y) 
and £2(2/) which parameterize the hadronic matrix elements of the weak currents between 
the states in the T doublet of bs system and those in the same doublet of cs system ap- 
proximately satisfy the normalization condition, £ x (y = 1) = 1 and £ 2 (y = 1) = 1, which 
is implied by the heavy quark flavor symmetry at the leading order of the heavy quark 
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expansion. 
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